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We discuss the propagation and fractionalization of localized charges on the edges of quantum Hall 
bars of variable widths, where interactions between the edges give rise to Luttinger liquid behavior 
with a non-trivial interaction parameter g. We focus in particular on the separation of an initial 
charge pulse into a sharply defined front charge and a broader tail. The front pulse describes an 
adiabatically dressed electron which carries a non-integer charge, which is times the electron 
charge. We discuss how the presence of this fractional charge can, in principle, be detected through 
measurements of the noise in the current created by tunneling of electrons into the system. The 
results are illustrated by numerical simulations of a simplified model of the Hall bar. 

PACS numbers: 71.10.Pm, 73.63.Nm 

I. INTRODUCTION 

Charge fractionalization, the appearance of quasiparticles which carry a fraction of the charge unit, is a remarkable 
effect that is found in certain quantum many-body systems with unusual properties. Well studied examples are 
the two-dimensional electron fluids of the quantum Hall effect, and ever since the basic theoretical understanding 
of the fractional effect was established it has been known that the fundamental quasiparticles in these systems are 
fractionally charged and satisfy fractional statisticfl^. Experimental studies of current fluctuations due to charge 
tunneling have indeed confirmed the presence of fractional charge carrierPH^. 

More recently it has been suggested that excitations with fractional charge may appear in one-dimensional systems 
described by Luttinger liquid theory^—. Here the fractionalization is linked to chiral separation of charges that are 
introduced in the systerrPifil, so that fractions of a unit charge move to the right and the left, respectively. In addition 
to theoretical works predicting this effect there have been suggestions of experiments that could detect the fractional 
chargeJiilti^, and an experiment has been performed that confirms the expected left-right asymmetry of the current 
associated with the injection of charges in the systerrfl^. 

There are however important differences between the fractionalization effect in the two systems, since the quantum 
Hall fluids are incompressible whereas a Luttinger liquid is gapless. This difference is of importance both for the 
question of uniqueness and of sharpness of the fractional charge. In the quantum Hall case, these properties of the 
quasiparticles follow from the topological properties of the fluids, while for the Luttinger liquid there seems not to be 
any unique value associated with the fractional charges, which are instead determined by the way these excitations 
are created. 

In a previous publication, three of us have examined questions concerning values and sharpness of fractional charges 
in Luttinger liquid^^ (see also Ref . [TH)) . The conclusion is that fractional charges in such systems can be sharp, not in 
an absolute sense, but in the sense that the charge fluctuations are indistinguishable from the background fluctuations 
of the ground state. These charges can take different values, depending on initial conditions, and we have in particular 
examined the difference between the situations where integer charges are introduced suddenly or adiabatically into 
the Luttinger liquid. 

The purpose of the present paper is to follow up this work by studying in some detail the time evolution of charges 
that are introduced as edge excitations in a quantum Hall bar, where the non-trivial Luttinger liquid behavior is due 
to interactions between charges on the two edges^*" We have performed explicit calculations of the time evolution 
of pulse shapes under transitions between regions with non-interacting and interacting edges, and focus in particular 
on what happens under (quasi-)adiabatic evolution from integer to fractional charges. 

Our Hall bar geometry is essentially the same as the one recently studied in Ref. [TBI where it was suggested that 
the charge fractionalization could be detected by noise measurements. Our analysis supports this claim, but our 
conclusions differ both on the expected values of the fractional charges, and on the optimal strategy for detecting 
them. 
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II. THE EFFECTIVE THEORY 



We consider a quantum Hall bar with a constriction at (Landau level) filling fraction ly — 1, see Fig. [T] The 
electrons are assumed to be fully polarized, so that the spin can be suppressed in the description. In the region of 
the constriction the two edges are sufficiently close to allow interaction across the sample, but sufficiently far apart to 
completely suppress charge tunneling between the edges. Outside the constriction, the separation between the edges 
is much larger, and interaction across the bar is completely suppressed. 
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FIG. 1: (Color online) Schematic illustration of the Hall bar with a constriction. At the point T an electron tunnels into 
the upper edge, shown in the form of sharply defined modulation of the edge (1). The pulse travels to the right, and when 
entering the constriction, the charge pulse is separated in a reflected and a transmitted pulse (2), with the transmitted pulse 
appearing in the form of a correlated modulation of both edges. This pulse is further split into a reflected one, moving within 
the constriction, and a transmitted one that propagates to the right outside the constriction (3). 



Electrons are assumed to be inserted by tunneling into the system at one of the edges in the region outside the 
constriction, and to travel with the edge current from the initial region into the constriction. The figure shows the 
qualitative picture for the propagation of the electron charge from the point of insertion. 

With no tunneling the charges on each edge are separately conserved and can be described by two independent 
fields p±, which measure the edge charge densities relative to their ground state values. The action can be separated 
in two parts. 



5* — 5*0 + Sint 

where Sq describes free spin-polarized electrons, and Sint is the electron-electron interaction, 

Sint = - ^ X! / dtdxidx2 [Px{xi,t) Vi{xi,X2) Px{x2,t) + Px{xi,t) V2{xi,X2) P-x{x2,t)] 

x=±-' 



(1) 



(2) 



where Vi is the interaction between charges on the same edge and V2 the interaction between charges on opposite 
sides. Symmetry between the two edges is assumed, and x is the linear coordinate in the direction of the symmetry 
axis. For a soft edge profile, and for sufficiently low energies, we may make the local approximation Va{xi,X2) = 
Va{xi)5(xi — X2), a = 1, 2. The resulting action is that of a Luttinger model 



S = ttH I dxdt ( dxQ dt^ — ^v{x) 



with variable parameters v{x) and g{x) given by 

v{x) — 



{d^<^f + g{x){d^Qf 



(3) 
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U+2k(Vl{x)~V2{x')) 
U+2kiVl{x)+V2{x))' 



V2{X) 



(4) 



where u is the Fermi velocity of the non-interacting theory. The two fields $ and 8 are related to the charge densities 
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through 



P±^\{d.^Td^Q) (5) 



The description in terms of the effective parameters v{x) and gix), as given by (4|, depends on the assumption of 
a smooth, quasi-adiabatic transition from the wide to the narrow part of the Hall bar. The expressions are then the 
same as with x-independent interactions Vi and V^^. A smooth modulation of the edge profile is desirable for the 
study of intrinsic physical properties of the system within the constriction, in which case effects that depend on the 
precise profile of the edges in the transition region are less important. From now on, we shall employ the effective 
model ([3]). We have checked numerically that for typical parameters, and pulse shapes, used in the below analysis, 
this effective model gives essentially the same result as the non-local microscopic theory defined by ([T]) and ([2]). The 
latter is presumably more accurate, but also computationally much more demanding. 

III. TRANSITION AND REFLECTION OF CHARGED PULSES 
A. Multiple reflection analysis 

It is convenient to change to new, g-dependent, charge density variables, defined by 

/± = ;^(5.$Tff5,e) (6) 
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They satisfy the following field equations. 



dt!±=TdAvf±)±\v^{f+ + n (7) 
2 g 

as can be derived from the action (|3|. In regions where g is constant f± define the two chiral components of the 
charge density, corresponding to the right and left moving parts of the edge fields. For the case 5 = 1 this separation 
in terms of the right and left moving components is identical to the separation of the total charge into components of 
the two edges, so that /± = p±. However, when g ^1 there is a difference between these two ways to decompose the 
total charge, and we have 

/± = ^(l±g)p+ + ^(lT5)P- (8) 

As a consequence there is for a purely right (left) moving mode (/zp = 0) a unique ratio between the charge densities 
of the two edges 

P±/Pt = ^, (/t = 0) (9) 

In a region with constant g a right(lcft) moving charge will thus decompose in two parts, with a charge distribution 
on the upper(lower) edge and a co-moving mirror image of the charge distribution on the lower(upper) edge, and with 
a fixed ratio {g + l)/{g — 1) between the two charges^S. 

The Luttinger model for the quantum Hall bar gives rigid constraints for charge transfer between two regions with 
different values of the interaction parameter g. This follows since, on one hand, the charge is separately conserved 
for each edge of the Hall bar, and on the other hand, the charge with a given chirality splits up in components on 
the two edges, which have a unique ratio determined by the value of g. As a result the transmission coefficient T and 
the reflection coefficient R for scattering of a charge on the boundary between two regions with different values of the 
interaction parameter are determined as 

T=^, R-'-^ (10) 
9 +9 9 + 9 

with g' as the value of the interaction parameter in the region of the incoming charge and g as the parameter in 
the region of the transmitted charge. In particular, for scattering from the non-interacting region [g' = 1) into the 
interacting region (g 7^ 1) we have T = 2g/(l -t- g) and i? = (1 — g) /(I -I- g), consistent with the results of Ref. One 
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FIG. 2: (Color online) Schematic illustration of charge reflections and transmissions in the discretized model. In a) a sharply 
defined charge pulse, represented by the thick black line enters a region D with a variable value of the interaction parameter g. 
On each step, where g changes, it is split into a transmitted and a reflected charge. The total transmitted and reflected charges 
can be separated in charge components, characterized by the number of internal reflections they have experienced during the 
transit through D. The transmitted component To, with no internal reflections, corresponds to the sharply defined front pulse 
of the transmitted charge. The refiected component _Ri, with one internal refiection, is broadened during the transit of D. 
It is clear from the illustration that the width of each charge component increases with the number of reflections. In b) the 
transmission and reflection at a single step is shown in a space-time diagram. The rules for reflection and transmission ( |10[ |, 
when applied at each step, gives rise in the continuum limit to the effective field equation of the system, as discussed in the 
text. 

It is instructive to consider a discretized model, as shown schematically in Fig.[2j for the transmission and reflection 
of charges in a region where g changes. The value of g (and of v) is then assumed to change in a stepwise fashion, 
with transmission and reflection of charge taking place at each step. The scattering at individual steps (as illustrated 
in Fig. [2]d)) is assumed to satisfy the relations (10 1 for the (local) transmission and reflection coefficients. 



The outgoing, transmitted and reflected charges can be decomposed into parts defined by a given number of 
reflections inside the interval where g changes, and an interesting point is that in the continuum limit of the discretrized 
function g{x), each of these charge components tends to a finite value, which takes a simple form when expressed in 
terms of the interaction parameter. With g = 1 in the region of the incoming pulse and g = exp(— 2r) in the region 
of the transmitted pulse, the expansions in number of reflections are expressed as 

2 24 720 8064 ^ 

i? = r--r3 + Ar5-ilr^ + — r9 + o(r") (ii) 

3 15 315 2835 ^ ^ ^ ' 
where the number of reflections corresponds to the power of F. The series sum up to T = and R = consistent 



with the asymptotic expressions given in (10 1. 

It is of interest to note that the effective theory described by the action ^ can be viewed as a local implementation 



of the relations ( 10 1 for transmission and reflection of charges. Thus we expect the continuum limit of the discrete 
model outlined above to give a faithful representation of the effective theory ([s]). To explicitly show this we consider 
the scattering of charge on a single step, as illustrated in Fig. ^p, in the form of a matrix equation, 

(12) 



where j± = ivf± define the charge currents. For an infinitesimal change in g, we have according to (10 1, T± = 
1 ± dg/2g and R± = ±dg/2g, and for corresponding infinitesimal changes in the coordinates x and t between the in 
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and out states, the scattering equations can be written as (see Fig.[2|D)), 

j+{x + dx,t + dt) = (1 + ^) j+(x,t) - ^ j_ (x + dx, 

j^ix,t + dt) - il-^)j^ix + dx,t) + ^U{x,t) (13) 
2,9 2g 

To first order in the differentials, and by use of the relation dx = vdt, the above equations can be re-written as the 
following differential equation 

dtj± = Tvd^j± + l,v — {j+- j-) (14) 
2 g 

and it is straightforward to check that it is equivalent to the field equation ([t]). 

An important point to note from the discretized model, clearly demonstrated in Fig. [2] is that the transmitted 
pulse has a distinct separation in a front pulse and a tail. Thus, with a sharply defined incoming pulse, as in Fig. [2] 
the leading transmitted pulse, which has no internal reflections, will be equally sharp. This is different from the other 
transmitted components, which are broadened by the reflections. With L as the width of the scattering region the 
contribution from n reflections has a width that is in fact proportional to nL, since each additional pair of backward- 
forward reflections will spread any part of the pulse over the width L, as is illustrated in Fig. [2j The net result is that 
the transmitted pulse obtains a long tail in addition to the sharply defined front pulse. This effect is present also in 
the continuum limit. We note that the transmitted front pulse has a shape that is insensitive to the functional form 
of g{x) and v{x), while the precise form of the tail will depend on these functions. 

The distinction between the front pulse and the tail is directly related to the discussion given in Ref. [TS] of how a 
charge that is adiabatically introduced into the system is separated into a local charge and a non-local charge that is 
evenly distributed over the system. The charge of the front pulse, Tg = is indeed identical to the charge found in 
Ref. [T3]for the local part of charge. The tail of the charge distribution, defined by multiple reflections in the transition 
region, will tend to the evenly distributed non-local charge in the adiabatic limit. For the reflected charge there is no 
component corresponding to the sharply defined front pulse, since all components are broadened by the reflections. 
For a smooth function g{x) the reflected pulse will therefore give a week, broad signal. 



B. Numerical simulations of pulse shapes 

Here and in the following we shall present results of numerical solutions of the field equation ([?]) with a discretized 
time coordinate, and with gaussian initial pulses at the point of tunneling. In this, and subsequent, simulations the 
shape of the edges are expressed in terms of error functions as 

with x the coordinate along the Hall bar and y the transverse coordinate, with W as the transverse width of the 
bar outside the constriction and Wc as the width within the constriction. The parameter b defines the longitudinal 
width of the transition region and i the length of the constriction. In dimensionless units the parameter values used 
in the evaluations are W = 10, Wc — 0.1 and b = 80. Below we shall consider two situations. In the first, where the 
tunneling occur at an external lead, we set i — 120, and in the second, where the tunneling is within the constriction, 
we set £ = 180. All other parameters are taken equal. Furthermore the width of the initial gaussian pulse of the 
tunneling charge is chosen as A = 0.6. A screened interaction between the charges on the two edges is assumed, and 
is here modelled by a Gaussian potential with damping length d « 0.8. The strength of the interaction is chosen to 
give g = 0.64 for the interaction parameter and Vc = 21 for the effective velocity within the constriction. 

All the above parameters are dimensionless and we need a suitable timescale for later evaluation of current and 
noise. Our choice is the propagation time between the endpoints of the constiction, i.e., the length of the constriction 
divided by Vc- The Luttinger parameter g{x) will not change until the distance between edges are of the same size as 
the interaction length. Thus the natural length of the constriction is not £, but an effective length £^g which we define 
to be the distance between the points where g{x) is halfway between 1 and and the value of g in the constriction. We 
will scale time with to = ics/vc and frequency with ujq = 2'!tVc/£cS in current and noise plots. 

In Fig. |3j we simulate the time evolution of a charged pulse incident on a single step in the width corresponding to 
a transition from a region with g = 1 at the far left to a region with g — 0.51 at the far right. The initial pulse is in 
the transition region split into a right moving, transmitted pulse and a left moving reflected one. The figure clearly 
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FIG. 3; (Color online) Numerical simulation of the scattering of an incoming pulse due to changes in the interaction parameter 
from g — 1 to g = 0.51. The figure displays the charge density at the upper (blue line) and the lower edge (red line) at different 
times. The lower (black) curve represents the changing value of the parameter g{x). Initially a pulse at the upper edge travels 
to the right in a region where g = 1. As it progagates through the transition region there are two distinct charges appearing 
on the lower edge. The mirror charge has the same shape as the pulse on the upper edge, but a smaller value and the opposite 
sign. The reflected part is much broader and its width depends on the length of the transition region. A close examination 
(see inset) reveals that the pulses traveling to the right divide into a sharp front pulse and a long tail arising from multiple 
reflections within the transition region. The shape of the tail and of the reflected pulse is determined by the profile of g{x) and 
v{x). The front pulse is independent of the shape of the functions, though its total charge is determined by the final value of 
g and its width on the final value of v. Since v decreases the front pulse gets compressed. 



displays the front pulse of the transmitted charge, vifhich is followed by a long, weak tail. Also the the broad reflected 
pulse is clearly visible. The two components of the front pulse, the main component on the upper edge and the mirror 
image on the lower edge, are both shown in the figure. 

We have applied two different mehods to estimate the charge of the front pulse from the numerics. The first method 
is to integrate the density distribution only in the region corresponding to the front pulse. Fig. 4(a) shows the sum of 
the density distributions at the upper and lower edge at the final time step of Fig. [3] We can identify the front pulse 
as the part of the density distribution where the density is positive, that is, for x > 25.6. A numerical integration 
in this region gives the value 0.7029 for the effective model and 0.7028 for the microscopic model, as compared to 
y/g = •\/0.5125 = 0.7159. It is interesting to compare these results to the total transmitted charge, which is given 
by the transmission coefficient, T — — 0.6777. The second method is to estimate the charge from a Gaussian 
fit to the front pulse. The blue line in Fig. |4(b)| shows the Fourier transform of the density distribution from the 
microscopic model in Fig. 4(a) We see that for A: — >■ the Fourier transform approaches the total transmitted charge, 
0.6777. The frontpulse is assumed to retain its shape through the transition region, we therefore fit a Gaussian to the 
Fourier transform for k > kco, where kco is a cut off. With the cut off equal to the width of the Fourier transform 
of the Gaussian we find that the charge is given by 0.7154. When varying the kco from 0.6 to 1.4 of the width of 
the Fourier transform of the Gaussian, the estimated charge varies from 0.7153 to 0.7162. This way of estimating the 
charge gives a better agreement with the predicted value of ^/g, than the first method. This is not surprising since. 
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with the parameter values we use, the step and pulse width are comparable in size. This means that the positive 
region is likely to have contributions from higher order reflections, and these will be included in the first estimate. 
In the second method only the high momentum components are used to fit the gaussian profile, and these are not 
expected to be sensitive to the finite width of the step. 
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FIG. 4: (Color online) Estimation of the charge of the front pulse, (a) Sum of the density distributions at the upper and lower 
edge at the final time step of Fig.|3] Only the rightmost part of the density is included, which corresponds to the rightmoving 
charge. The inset helps to distinguish between the front pulse and the tail, (b) Fourier transform of the density distribution 
of the microscopic model (blue curve), shown with a Gaussian fit to the curve (in green). The red rings are points used to 
determine the fit parmeters (only every 20th point is shown). From the inset we clearly see that the blue line approaches the 
value of the total transmitted charge for fc — > 0, while the green line approaches y/g. 



IV. CURRENT FLUCTUATIONS AND FRACTIONAL CHARGES 



By measuring current fluctuations, one can gain information about the presence of fractional charges in the system. 
Ref. ITT] has considered scattering of charges on an impurity in a system with an effective Luttinger parameter g, 
where noise is produced in the backscattered current. The expected value for the scattered charge, derived from 
the a; — >■ limit of the Fourier transformed noise, was found to be g times the electron charge. In Ref. [13l which 
considered the Hall bar geometry discussed above, it was suggested that the fractional charge (1 — .g)e/(l + g), given 
by the reflection coefficient in (10 1, could be observed by measuring the tunneling noise in the reflected current on 
the lower edge (see Fig. [T]). However, secondary pulses created from reflections on the second boundary will in this 
case neutralize the charge to form an extended reflected signal with zero total charge, and the information about the 
presence of fractional charge therefore has to be extracted from the w 7^ part of the noise. 

We too consider the current noise due to tunneling, but focus on the above distinction between the front pulse, 
and the broadened tail, of a charge transmitted through the constriction. The front, which is the robust part of 
the pulse, carries the charge ^/ge, and we take this as the natural definition of a fractional charge value within the 
constriction. The expressions for the noise are found by a simple, semiclassical approach, and they essentially agree 
with the expressions given in Ref. [T31 Although it might be more difficult to measure, we shall consider the noise also 
in the current inside the constriction, since it provides a clear signature for fractional charge. In the last section we 
shall propose an alternative tunneling geometry where there is a clear signature in a more readily measurable current. 



A. Noise in tunneling currents 

The tunneling is described as a Poisson distributed sequence of events, each characterized by a sharply defined 
current pulse Is{t), with J Is{t)dt = e. The total current at time i, restricted to a (large) time interval T, at a 
point close the to the point of tunneling, is a distribution of N pulses superimposed on the background edge current. 
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hot = hun + lo, with 

N 



n=l 

where i„ denotes the instant of a tunnehng. The tunneUng events are assumed to occur randomly within the given 
time interval, thus giving a constant average current 

dtnlsit - tn) ^ e- . (17) 



{Itu7l) — / 
n=l ^ 



This average is held fixed, and in the following the limit T, A'^ — > oo is taken whenever convenient. 

The noise of the current is defined as the Fourier transform of the current-current correlation function 

/oo 
dte-*[(/W/(0)>-(/(t)>(/(0)>]+c.c.. (18) 
-oo 

With no correlation between the tunneling and background currents we have S{io) — Stun{^) + 'S'o(w), and we focus 
primarily on the fluctuations in the tunneling current. The assumption of randomly distributed tunneling events 
implies that the expectation values can be found by independently integrating over the tunneling times i„. 



Stuni'^) = hm 

T— i-oo 



1 I'ca I'l 
I r°a rT rT 



POO 



dte'^' {Itunf 



+ c.c. 



= 2^i,{oj)h{~uj) = 2e|/,(t^)|2 (19) 

where /^(a;) is the Fourier transform of Is{t) and /s(w) is the normalized profile function, defined by Is{uj) — eja{u)). 
This definition gives /s(0) — 1 and therefore Stund^ — >■ 0) = 2e consistent with the expected noise-current 

relation in the non- interacting region (g = 1). As shown by the above expression, the noise can be viewed as a single- 
pulse eff'ect, since it is determined by the Fourier transform Igiuj) of the charge pulses generated by the tunneling. 



B. Noise in the reflected current 



Due to the linear propagation of the charged pulses through the Hall bar, the current fiuctuations at other points 
along the bar are related in a simple way to the fluctuations at the initial point. Note that no additional noise is 
produced by the constriction itself due to the adiabatic transition from the outside to the inside. We consider flrst 
the fluctuations in the reflected current on the lower edge (Fig. [T]). The relation between the current at this point 
and at the initial point can be written as 

/oo 
dt'GR{t-t')Itun{t') (20) 
-OO 

with Gji{t) as the propagator between the two points. The Fourier transform takes the form 

iRiuj) = R[uj)itun{i0) (21) 

where R{io) — J dtexp{iLdt)Gu{t) is the frequency dependent reflection coefficient, with i?(0) as the full reflection 
coefficient. It follows from the definition of S{u!) that the noise at the point where the reflected current is measured 
is given by 



Sb.{^) = |-R(w)pS'f„„(a;) 
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= 2e|i?Hn/,H|2 (/,„„) (22) 

The mean reflected current is {In) = R{0) {Itun) ■ However, the transmission coefficient for the Hall bar across the 
constriction is T = 1, since we have the same value g = 1 on both side^. This means that the reflection coefficient is 
R = R{0) — and therefore there is no contribution to the average edge current from the tunneling, (//j) = 0. Also 
the tunneling contribution to the noise vanishes in the uj limit. 
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FIG. 5: (Color online) The figure displays the result from a numerical evaluation of the tunneling contribution to the noise, 
Sr{uj), of the reflected current on the lower left edge of the Hall bar. The value of the interaction parameter is g = 0.64 and the 
frequency scale ojo = 2TTVc/£ef! is given by the effective velocity Vc = 21 and length i^g = 77 of the constriction. The green curve 
is the contribution from the first refiected pulse only, which is reflected at the left boundary of the constriction. The numerically 
determined value for the — > limit of this contribution fits with high precision the theoretical value [{g— !)/(<?+ 1)]^ ~ 0.048. 
The decay length of the curve is inversely proportional to the width of the reflected pulse, and thereby to the width of the 
transition region where g{x) changes. The blue curve represents the full noise function, where contributions from the secondary 
pulses are included. These are pulses that arise from multiple reflections inside the constriction. The value that is found for the 
full noise function in the limit cj — >■ is to high precission, consistent with the expectation that charges of the reflected pulses 
add up to zero. The oscillations in the blue curve are due to the time shift between pulses reflected at the two boundaries, and 
the regular form of the curve is due to the symmetric form of the constriction. 

The vanishing of the average current can be viewed as due to cancellations between charge components associated 
with repeated reflections between the boundaries of the constriction. In Ref. [13] it has been suggested that one 
can nevertheless extract information from the noise, by exploiting the time delay of the secondary reflected pulse, 
either by filtering out pulses from the second boundary or by measuring the noise function S'_r(w) for non- vanishing 
frequencies uj. The noise function Sii{uj) is here obtained by a numerical Fourier transform of the time dependent, 
reflected current, and is measured relative to 2e (Itun)- The result is displayed in Fig. [5| which shows both the full 
function, and the function defined by including only the contributions from the first reflected pulse. As expected, the 
full function tends to when w — > 0, and it shows an oscillatory behavior due the interference between contributions 
from repeated reflections. The decay length of the curve is inversely proportional to the time width of the flrst 
reflected pulse, which is in turn proportional to the width of the transition region where g changes. Even though 
the oscillatory function contains some information about the contribution from the first reflected pulse, shown as the 
smooth green curve in the figure, and thereby about the charge carried by this part of the refiected signal, to extract 
a good numerical value from a corresponding experimental curve may seem difficult. Note also that the current and 
noise shown in the plot in reality sits on the top of the background edge current, which has both non- vanishing average 
value and noise. 



C. Current noise within the constriction 

As pointed out, measurements on the reflected current is hampered by the signal being weak, and broadened by 
the width of the transition regions where g changes. In fact, in the adiabatic limit the reflected pulses completely 




■ First pulse 

■ Full current 
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disappear into background current. For the transmitted current, where the front of the pulse is not broadened in the 
transition region, the situation is quahtatively different. We consider therefore next the current and current noise 
within the constriction, assuming simply that these can also be subject to measurements. Since the total current, 
rather than the currents on the separate edges will be more accessible to measurements, we will focus on this quantity, 
see Fig. [6] 
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FIG. 6; (Color online) Time dependence of the current due to a single tunneling event, numerically evaluated at a point within 
the constriction. The time scale to — iad/vc is given by the effective velocity Vc = 21 and length ics = 77 of the constriction. 
The front pulse is seen as the sharply defined peak, followed by the tail which is so small that it is barely visible. The smaller 
pulse at f ~ 1.8 carries the charge reflected at the right end of the system (pulse 3 in Figure 1). The charge of the pulse is 
negative, but it gives a positive contribution to the current since it is propagating to the left. It is considerably broader than 
the first pulse as a consequence of the reflection. At later times more pulses will pass, but their values decrease and their widths 
increase with each reflection, making them rapidly less significant. There is also a small second refiected pulse, at t ^ 2.8, that 
is not visible in the plot. 

The current in the constriction is related to the initial tunneling current in precisely the same way as the reflected 
current. We write it as 

/oo 
dt'GT{t-t')Itun{t') (23) 
-OO 

with G*t(^ — i') as the propagator from the initial point to the point where the current is measured. We similarly 
have for the Fourier transform 

/oo 
dte^^'Grit), (24) 
-OO 

and for the current noise 

Sriio) = |T(^)p^t™(^) = 2e|r(c^)|2|/,(t^)|2 . (25) 

Since the transmission coefficient is unity, T = T{0) = 1, the average current is simply (It) = {hun)- Therefore a 
naive application of the ratio between the noise and average current in the w — > limit would give e as the charge, 
rather than a fractional value. 

However, in the same way as for the reflected current, this trivial result is caused by multiple reflections between 
the boundaries of the constriction. To compensate for the effect of these multiple reflections, which mask the presence 
of non-integer charges in the limit a; — >■ 0, it is natural also here to consider the noise for non-vanishing lu. The 
broadening of the signal, and the reflections between the two boundaries, will affect the low frequency part of the 
noise, while the form of the front pulse will shape the higher frequency part. Thus the transmission coefficient |T(w)|, 
which tends to 1 for small w, is expected to approach for large oj. The effect on the noise is shown in Fig. [71 where 
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FIG. 7: (Color online) The numerically evaluated noise function Sri^jj) at a point within the constriction, plotted as a function 
of frequency to. The plot to the right shows details of the left plot for small to. The green curve represents the result when only 
the first pulse traveling through the constriction is included, while the blue curve gives the result for the full current. At high 
frequencies there is overlap between the green and the blue curves, indicating that the high frequency part of the noise function 
is governed by the front pulse alone. The decay length of the noise function is then determined by the width of the front pulse. 
For decreasing values of u both curves approach the value g = 0.64 (indicated by the red, dashed line), corresponding to the 
square of the front pulse charge, but then deviate from this for even smaller values of lj. For the green curve the deviation can 
be explained as due to charge contribution from the long tail of the front pulse. For the blue curve the oscillations are caused 
by the charge reflections between the boundaries of the constriction. The damping length of the oscillations is determined by 
the width of the transition region where g{x) changes. The frequency scale is given by the effective velocity Vc — 21 and the 
effective length of the constriction /eft = 77 through ujq = 27r«c//cft . 



the broadening and multiple reflections give rise to oscillations for small w, while the profile of the noise function for 
larger u! is determined by the front pulse. 

In this case the time extension of the front pulse determines the (inverse) width of the noise function, which is 
therefore much wider than the noise function of the reflected current, and provided the relevant time scales are well 
separated, the value of the fractional charge can then be extracted for values of uj where the oscillations are strongly 
damped. Thus, in this case the value of the fractional charge y/g is (in principle) easier to extract since the signal is 
much clearer. Note, however that 5t(w) depends quadratically on the transmission coefficient T(uj). The peak value 
of ST{uj)/{2e {I tun)) is therefore closer to the quadratic value g than to the charge value ^ of the front pulse. 



V. FINITE TEMPERATURE EFFECTS AND BACKGROUND FLUCTUATIONS 



In the expressions for the noise used so far, we have implicitly assumed zero temperature, since only tunneling 
into the Hall bar is assumed. The effect of finite temperature can be taken into account by assuming that both 
tunneling into the Hall bar (creation of a charge) and tunneling out of the Hall bar (creation of a hole) can take 
place, with the relative probability of these two types of events being determined by the Gibbs factor exp(/3A/x), 
where A/i — ^2 — fJ-i = sVtun is the difference in electrochemical potential between the edge of the Hall bar and 
the tunneling reservoir respectively. To be more specific, if we assume the temperature to be high enough for the 
Boltzmann distribution to be valid, the number of electrons, iV+, and holes, A^-, that are randomly injected into the 
system are given by 

^ ^ exp(/3^2) ^ N exp(i/3A/i) 

exp(/3/^i) + exp(/3^2) 2 cosh(i/3A/i) 
^_ ^ exp(/3/.i) ^ N_ cxp(-^/3Am) 

exp(;3^i) +exp(/3^2) 2 cosh(i/3A^) ' 
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corresponding to the mean current 

fiV . - N ) N 1 
(Itun) = e } +^ > = e- tanh(-/3AM) . (27) 

If for simplicity we assume the profile functions for charges and holes to be the same, the expression for 5't„„(cij) will 
however remain unchanged, giving a temperature dependent relation between the noise and the mean value of the 
tunneling current, 

Stuni^) = 2e|/,(c^)p coth(^^A^) (/i„„) . (28) 



Since the system is linear, the same factor coth(i/3A/i) will modify the noise in the reflected current (Eq.(22)) and in 



the transmitted current (Eq.(25l). 

The total noise in the current is then a sum of the background noise, without the tunneling, and the noise in the 
tunneling current. The background noise is most easily determined in the bosonized theory. Assuming g{x) to be 
sufficiently smooth, and the velocities sufficiently small, for the pulses to be subject to an adiabatic variation in the 
value of g, the fluctuations are the same as for fixed g. In this case we may use the following momentum expansion 
of the current, in terms of the bosonic creation and annihilation operators bq and 




/(x, t) = I J ^ sgn(g) ioq [6,e'(«--'*) + fet e-^?--,*)] . (29) 



Lg_ 
2TT\q\ 

The system is here assumed to be confined to a ring of length L with g as a constant, and the relation between 
frequency ujq and momentum q is ujq = v\q\, with v as the velocity of the edge current. The expression for the noise 
is then, 

.2 ,,,2 



Soico) = /rfie-*^^^(6,6t+6t5^)cos(a.,t) 



9 

o2 



du)q LOq C0th( - /Sw^ ) [(5(lxl — LOq) + 5{uj + Uq)] 



f , , /I 

= / dWg a;„ coth(- 

Jo 2 

= ge"- coih{\pLo) , (30) 

where we have taken the continuum limit L oo and assumed a finite temperature Bose-Einstein distribution 
{b^jbq'j = (exp(/3wq) — 1)""'^. It is interesting to note that the background noise, which is quadratic in the electron 
charge, is linearly renormalized by the interaction parameter g. This means the variation with g is consistent with 
the picture of an adiabatic change of the electron charge e — ^ y^e, from the integer to the non-integer value, when 
the interaction is turned on. 



VI. TUNNELING WITHIN THE CONSTRICTION 



As previously discussed, the fractional charge y/ge, associated with adiabatically dressed electrons within the 
constriction, can in principle be detected in a measurement of the current noise in the constriction. Outside the 
constriction there is an indirect and less clear relation between the charge fractionalization and the noise in the 
reflected current as well as in the transmitted current. In the transmitted current, outside the constriction, there 
is in fact no information about the fractional charge in the sharply defined front pulse, since the strength of this is 
redressed to the value 1 when g{x) regains this value. Therefore the information lies, as is the case for the reflected 
current, in the form of the low frequency part of the noise, which depends on how the repeated reflections inside the 
system shape the signal. 

This motivates us to consider a different tunneling scenario. Instead of coupling to an external reservoir, we assume 
the distance between the edges of the Hall bar in the constriction to be sufficiently small to have a non-vanishing 
probability for charge to tunnel between them. A potential difference A/i between the two edges will then introduce 
a tunneling current. As a basic assumption we take that each tunneling event corresponds to moving one electron 
between the edges. We use the same picture of the tunneling current as before, now with denoting the number 
of transitions of one electron from the lower to the upper edge in a given time interval T and as the number 
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of transitions in the opposite direction. With these assumptions the average tunnehng current, at non- vanishing 
temperature, is as before given by Eq.(27). Charge conservation on each edge implies that this is also the average 



current outside the constriction to the right, while the current to the left is the same in absolute value, but with 
opposite sign. 



N 1 

{Itun) = {Iright) = - {Ileft) = 6— tanh(-/3A/x) 



(31) 



The expressions for the noise are essentially the same as before, with the noise in the tunneling current given by 
Eq.(l28|, 



1 



Stuniio) = 2e\fs{uj)\^coth{-l3Afi) 



(32) 



with fs(to) again denoting the profile function of the tunneling charge pulse. The noise of the right moving current 
outside the constriction we write as 



Sright{i^) = 2e|THght(a;)|^|/5(a;)|^coth(^^A/i) [hun) 



(33) 



with Trighti^) as the Fourier transform of the propagator from the point of tunneling to a point outside the constriction 
to the right, see Fig. |8] The noise in the left moving current is the same. 
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FIG. 8: (Color online) Numerical results for the zero temperature noise outside the constriction, reproduced from simulation 
of the charge propagation from one tunneling event within the constriction. The green curve also here includes only the 
contribution from the first pulse, and the blue curve is the noise produced by the full current, i.e., with secondary pulses 
included. The noise function is very similar to the noise function in FigjT] for the case of tunneling taking place on the outside 
of the constriction. Again, the noise of the front pulse, for small lu, dips below the value g = 0.64 (dashed red line), due to 
contributions from the long tail, and the noise of the full current shows for small w oscillations due to the efi'ects of secondary 
pulses. The frequency scale is given by the effective velocity Vc ~ 21 and the effective length of the constriction i^g = 117 

through Wo = 2TVVc/ieB. 



To get a qualitative understanding of how these expressions relate to fractional charges within the constriction, 
we first note that the sudden transition of an electron from the lower to the upper edge will create a right moving 
pulse within the constriction of strength g and a left moving pulse of strength —g. These values are again determined 
by charge conservation on each edge combined with the fixed ratio between the charges on the two edges for each 
chiral component. The transmission of this pulse to the outside gives rise to the right moving front pulse with charge 
determined by multiplication with the inverse dressing factor 1/y^. Thus the charge is ^/ge, precisely the same 
charge as for an electron that moves from the non-interacting region into the constriction. The total charge of the 
right moving pulse is however identical to the electron charge e, due to charge conservation on the upper edge. This 
charge includes the contribution from the tail and the secondary pulses created by reflections within the constriction. 

The expected form for the transmission coefficient Tright{^) is then much the same as for Tiu), as previously 
discussed. For small lo the form is determined by the broadening and multiple reflections within the constriction. 
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Again the limit w — gives Tright{0) — 1, consistent with the mean value of the current not depending of g, 
(Iright) — (Itun) ■ For larger values of u the transmission coefficient approaches the limit ^Jg determined by the front 
charge. The result is that the noise shows essentially the same frequency dependence here as in the previous case when 
measured within the constriction. The main difference is that in this case, with tunneling inside the constriction, 
the information about the fractional charge can be extracted by current measurements on the upper edge outside the 
constriction, where the measurement can more easily be performed. 



VII. SUMMARY AND OUTLOOK 



In this paper we have shown how an effective Luttinger model, Eq. (|3]), can be used to simulate the low-energy 
dynamics of charged edge pulses on a Hall bar with a smoothly varying width. In particular, we have studied 
the reflection and transmission of such pulses from smooth contractions, described by a space dependent Luttinger 
parameter g{x), and also analyzed the current noise due to tunneling processes. Our theoretical analysis shows that a 
charge e pulse transmitted from a non-interacting region (g = 1) into an interacting region characterized by g ^ 1, is 
composed of a sharp, charge e-^g, front pulse that suffers no backscattering, and a tail, with a width depending on the 
velocity and the size of the transition region, due to repeated reflections. The reflected pulse, on the contrary, is not 
expected to show any sharp feature. All these theoretical predictions were conflrmed by numerical simulations which 
clearly show both the presence of a well deflned e^/g transmitted pulse with a broad tail, and a broad reflected pulse. 
We have also analyzed the current noise due to tunneling both from an external reservoir, and within the constriction, 
using a simple formalism based on the bosonized effective action (|3|) . For the former case we concluded that the e^/g 
charge of the front pulse can be extracted from the noise function A'(w) provided the latter can be measured at finite 
uj inside the constriction. If the tunneling takes place within the constriction, the fractional charge can be extracted 
from measurements of S{uj) outside the constriction. 

A line of investigation made possible by this work is to use the methods developed in Ref. [T^to study the properties 
of quantum noise due to the propagation of fractionally charged pulses in regions of varying g. The most interesting 
implication of this work, however, is the possibility to detect fractional charge by inducing electron tunneling 
within a constriction, and measuring the current noise on an external lead. Both the geometry, and the predicted 
value for the charge differ from earlier proposals. In this paper we have concentrated on the conceptual aspects of 
the problem, but it is clearly of great interest to examine if experimental realizations are possible. This may involve 
more detailed numerical simulations with realistic experimental configurations. In particular, it is important to find 
out if the effective Luttinger model can describe a realistic experiment. 
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